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I . INTRODUCTION 


Blast  loading  phenomena  are  major  contributors  to  both  engineering 
and  military  operational  environments.  Scientific  interest  was  evident 
as  early  as  the  nineteenth  century! >2,  continuing  effort  to 

characterize  these  phenomenological  effects  on  structures,  hiomans,  and 
military  targets  includes  the  work  of  Sperrazza^,  and  0.  T.  Johnson‘S, 
and  Baker!.  As  a continuing  part  of  this  effort  the  current  authors 
propose  a model  for  evaluation  of  blast-target  interaction.  The  work 
of  several  predecessors  at  the  Ballistic  Research  Laboratories  has  been 
influential  in  forming  the  proposed  model.  Consequently,  the  present 
report  forms  a kind  of  status  report  on  blast  vulnerability  modeling 
within  the  US  Army's  Ballistic  Research  Laboratories. 

The  first  concept  contributing  to  the  current  work  is  an  iso-damage 
model  using  a target,  a damage  level,  a blast  pressure,  and  a total  blast 
impulse  to  characterize  blast  effects^. 

The  second  contributing  analysis  is  the  work  of  Youngdahl^  in  which 
it  was  found  that  for  rigid-plastic  structures,  subjected  to  transient 
loading,  the  extent  of  damage  could  be  characterized  by  a two  parameter 
set,  i.e.,  those  of  total  impulse  and  a characteristic  time.  Both  are 
integrals  of  load  data.  This  was  an  important  achievement  and  an  indi- 
cation that  similar  methods  should  yield  similar  results  for  the  blast 
iso-damage  problem. 

An  important  third  contribution  to  the  authors'  thinking  was  the 
blast  damage  scaling  of  0.  T.  Johnson^,  That  approach  proposed  that 
level  of  damage,  target  and  explosive  source  should  combine  to  define  a 
unique  scaling  law  for  blast  damage. 

In  this  work  the  objectives  of  the  authors  were  to:  (1)  minimize 

subjectivity  of  blast/daraage  relationships;  (2)  to  simplify  data 
reduction;  and  (3)  to  improve  accuracy  of  blast  damage  analysis.  The 


^yilfred  E.  Bakev,  Explosions  in  Air,  Univ.  of  Texas  Press,  Austin, 
Texas,  1973. 

2 

F.  B.  Porzel,  "Introduction  To  A Unified  Theory  of  Explosions  (UTE) , " 
NOL  TR-72-209 , 14  September  1972,  Unclassified. 

2 

J.  Sperrazza,  "Dependence  of  External  Blast  Damage  to  the  A-25  Air- 
craft on  Peak  Pressure  and  Impulse,"  Ballistic  Research  Laboratories 
Memorandum  Report  575,  September  1951,  AD  # 378275. 

4 

0.  T.  Johnson,  "A  Blast-Damage  Relationship,"  Ballistic  Research 
Laboratories  Report  No.  1389,  September  1967,  AD  it  388909. 

5 

C.  K.  Youngdahl,  "Correlation  Parameters  for  Eliminating  the  Effect 
of  Pulse  Shape  on  Dynamic  Plastic  Deformation."  Journal  of  Applied 
Mechanics , pp.  744-752,  September  1970. 
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roles  of  the  prior  works  cited  have  been  to  provide  a basic  model 
(pressure  and  impulse*),  to  provide  an  additional  analysis  tool  (impulse 
and  a kind  of  moment  integral),  and,  to  provide  a means  of  validation 
of  the  process  (the  0.  T.  Johnson  scaling  law).  In  presenting  this 
work  the  proposed  model  is  first  defined  and  then  related  to  previous 
models.  Subsequently,  it  is  applied  (along  with  the  traditional  P-I 
model)  to  niiraerical  simulations  of  the  blast/damage  process.  The 
authors  believe  that  the  results  justify  wide  application  of  the  current 
mode 1 . 


II.  DESCRIPTION  OF  THE  MODEL 
The  proposed  model  is  this ; 

1.  Total  impulse  per  unit  area  delivered  by  the  blast  is  defined 
t = t^ 

'total  '/ 

t = 0 


where  time  start's  at  load  arrival  and  stops  at  the  end  of  positive 
pressure  phase,  (t^)  while  P(t)  is  the  pressure  delivered  at  the 
target. 


2.  Two  numbers,  representing  idealized  static  and  impulsive 

load  asymptotes,  which  cause  a specified  level  of  damage,  are  defined  as 
the  critical  pressure  (P^^.)  and  critical  total  impulse 

3.  A modified  pressure  is  defined  by  the  integral  expression 


P ^ 


P(t)dt 


fr 


stop 


start 


(2) 


where  t ^ * is  the  time  that  P(t)  first  exceeds  P and  t ^ 

start  . ^ cr  stop 

is  the  time  that  P(t)  drops  below  P^^  (see  Figure  1), 

4.  The  fundamental  blast/damage  model  is  expressed  in  the  form 


^Impulse  Zn  this  report  is  impulse  per  unit  area  unless  stated  otherwise. 
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Figure  1.  Critical  Pressure  and  Time 
Criteria  of  Proposed  Model. 
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CP  - P )(I^  ^ T - I ) = CONSTANT  C3) 

^ cr  total  cr-" 

5,  The  proposed  model  may  be  used  when  there  is  no  target  data 

(i.e.,  when  P and  I are  not  known)  to  evaluate  relative  blast/ 

^ * cr  cr  ^ 

damage  potential  by  writing  equation  (3),  with  = 1^^  = 0,  as 

P • I ^ = DN  (4) 

total 


where  DN  stands  for  "damage  number,"  a characteristic  of  the  blast 
pressure  only. 


Discussion  of  the  model  in  this  report  is  directed  toward  those 
levels  of  sophistication  which  the  authors  believe  justified  by  the 
past  and  current  practices  of  the  military  community.  There  is  no 
effort  to  place  the  presentation  in  the  context  found  in  the  physical 
science  studies  of  blast  phenomenology. 


III.  THE  MODEL  AND  PRIOR  WORK 

The  work  of  Sperrazza^  may  be  characterized  by  definition  of  an 
iso-damage  criteria  (see  Figure  2) 


P ) (I 
cr"^  ^ total 


- 'cr> 


CONSTANT 


where  Pp  is  peak  pressure  and  the  subscript  cr  identifies  target 

characteristics  of  damage.  Equations  (3)  and  (5)  are  of  exactly  the 
same  form  and  only  the  peak  pressure  of  equation  (5)  differs  from  the 
modified  pressure  of  equation  (3).  Equation  (5),  i.e,,  the  historical 
method,  carries  an  implicit  assumption  that  the  pressure  time  history 
P(t)  has  the  form 


P(t)  = Pp0(t,x,y,z,  ...) 


(6) 


hence  for  all  conditions  of  nature  0 has  little  influence  in  characterizing 
failure.  That  implicit  assumption  has  been  very  successful  in  the  hands 
of  experienced  blast  vulnerability  analysts.  It  has  also  given  rise  to 
some  difficulty  in  relating  test  data  to  ideal  theoretical  cases  and 
has  made  comparison  of  significantly  different  types  of  blast  loading 
difficult  (e.g.  Pentolite  versus  Fuel-Air  Explosive) . it  has  presented 
a source  of  error  in  data  reduction  due  to  the  need  for  subjective 


selection  of  peak  pressure 
that  the  work  of  Youngdahl 


§ 


It  is  in  removing  this  technical  problem 
provides  guidance. 


Youngdahl  proposed  the  use  of  a two  parameter  set  to  characterize 
the  damage  process  for  a rigid  perfectly-plastic  target  subjected  to 
transient  loading.  Like  Sperrazza  he  uses  total  impulse  as  one  para- 
meter. But,  instead  of  peak  pressure  he  uses  a characteristic  time, 
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Pressure 


Impulse  Per  Unit  Area 


Figure  2.  Pressure-Impulse  Iso-Damage  Model. 
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t , defined  by 


- A 
t = 


(t  - t^^)P(t)dt 


P(t)dt 


(7) 


where  t is  the  time  at  which  damage  starts  and  t is  the  time  at 
St  * sp 

which  damage  stops  (see  Figure  3) . Note  that  equation  (7)  is  a 

normalized  first  moment  of  the  loading.  Given  that  the  procedures  for 

blast  damage  analysis  in  the  military  community  already  provide  for  use 

of  the  form  of  equation  (5) , it  is  desirable  to  convert  equation  (7)  into 

a modified  pressure  rather  than  a time.  This  was  done  by  defining 


P = I/t 


(8) 


where 


I 


A 


(9) 


The  result  of  substituting  equations  (7)  and  (9)  into  equation  (8) 
is  similar  in  form  to  that  stated  in  equation  (2)  where  t^^  is 

approximated  by  and  t^^  is  approximated  by  The 

consequence_of  equation  (2)  (and  equation  (8))  is  that  the  number  pair 

I , and  P is  obtained  without  subjective  judgement.  Caution  must 
total 

be  exercised  when  attempting  to  analyze  data  when  the  function  P(t) 
repeatedly  changes  sign,  because  equation  (1)  implies  the  P-I  relation- 
ship applies  only  for  pressure  pulses  during  the  first  positive  over- 
pressure phase.  In  general,  multiple  zero  pressure  crossings  of  the 
loading  function  would  indicate  a repetitive  loading  or  quasi-steady 
state  forced  vibration  problem  rather  than  a pulse  loading  situation 
to  which  the  P-I  technique  should  be  limited. 


IV.  DISCUSSION  OF  THE  DAMAGE  NUMBER  MODEL 

Equation  (4)  is  the  form  in  which  the  proposed  model  must  be  used 

when  P = I =0.  This  also  implies  that  t . = 0 and  t = 

cr  cr  start  stop 
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t^,  i.e.,  that  all  of  the  positive  pressure  history  is  used  to  compute 

F,  because  the  starting  and  stopping  times  for  the  integrals  used  to 
compute  F are  associated  with  the  damage  pressure  threshold  (P  ) . We 

4 ^ 

must  now  refer  to  the  work  of  0,  T„  Johnson  to  establish  the  relative 
validity  of  equation  (4). 

0.  T.  Johnson's  work  relates  a normalized  nondimens ional  distance  at 
which  a given  level  of  damage  can  occur  (represented  in  Figure  4 by  C^, 
a factor  inversely  proportional  to  range)  to  the  quantity  of  explosive 
used.  As  is  clear  from  Figure  4,  many  datum  points  confirm  a scaling 
law  with  a slope  in  the  log-log  plot  equal  to  0.435.  In  Figure  5 a com- 
parison is  made  between  Johnson's  predictions  and  predictions  based  on 
the  Ballistic  Research  Laboratory  Pentolite  data®  using  equation  (4). 

The  predictions  are  based  on  a reference  point  for  some  level  of  damage 
(triangle  points  in  Figure  5)  having  known  characteristics  of  explosive 
charge  mass  and  standoff  distance  (range).  For  the  four  reference  points 
in  Figure  5 predictions  were  made  for  the  range  associated  with  a charge 
mass  of  100  kg  and  slopes  (in  the  log-log  plot)  for  the  DN  = constant 
predictions  determined.  The  mean  value  of  the  DN  slopes  obtained  in 
Figure  5 is  0.442  which  is  within  one-third  of  one  standard  deviation  of 
that  obtained  by  Johnson.  The  author's  interpret  the  agreement  between 
Johnson's  prediction  and  the  DN  prediction  to  be  a form  of  validation  of 
equation  (4),  the  least  sophisticated  form  of  the  proposed  model,  with 
the  Ballistic  Research  Laboratory  experimental  data  base. 

The  way  in  which  the  side-on  pressure  and  reflected  pressure  bracket 
Johnson's  data  (see  100  kg  in  Figure  5)  gives  an  indication  of  why  sub- 
jective estimates  of  damage  can  be  insensitive  to  details  of  the  damage 
mode.  Clearly  the  scaling  laws  for  both  side-on  and  reflected  pressure 
have  very  similar  forms,  with  only  the  initial  scale  requiring  adjustment. 


V.  STATIC  LOADS  AND  IDEAL  IMPULSIVE  LOADS 

The  idealized  limit  cases  in  which  engineering  load  analysis  are 
imbedded  are  the  static  load  and  the  pure  impulse  (P(t)  ->  “, 

duration  0).  There  is  no  conceptual  difficulty  in  applying  equations 
(1)  and  (2)  to  static  loading,  but  for  purely  impulsive  loading  some 
doubt  arises.  If  the  positive  duration  approaches  zero,  then  P may  not 
be  defined.  A careful  examination  of  the  application  of  equations  (1) 
and  (2)  to  a single  degree  of  freedom  oscillator  excited  by  a Dirac  delta 
function  shows  that  the  value  of  impulse  captured  is  defined  but  a 
value  for  F is  infinite.  This  is  the  expected  result  and  leads  to  the 
conclusion  that  the  method  has  generality.  But  in  order  to  ensure 


H.  J.  Goodman^  "Compiled  Free-Air  Blast  Data  on  Bare  Spherical 
Pentolite,"  Ballistic  Research  Laboratories  Report  No.  1092, 
February  1960,  AD  # 235278. 
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Figure  5.  Free-Air  Charge  Weight  Vs.  Range  for  Various  Damage  Numbers. 


greater  applicability,  the  authors  have  come  to  the  conclusion  tha^ 
the  use  of  the  model  should  probably  be  restricted  to  cases  where  P 
can  be  calculated  directly  without  access  to  sophisticated  mathematics. 
As  a practical  matter,  the  model  is  proposed  only  for  those  cases 
where  it  is  possible  to  determine  a critical  pressure  such  that  a static 
pressure  of  magnitude  causes  damage  equal  to  that  of  interest  in 

the  blast  loaded  case.  In  such  cases  P can  always  be  determined. 


VI.  CALCULATIONS  TO  TEST  THE  MODEL 

The  concept  of  the  model  is  of  value  if  it  achieves  any  one  of 
three  objectives:  fl)  that  it  allows  the  removal  of  subjective  judgement 

from  blast-damage  analysis  and  testing;  (2)  that  it  improves  data  re- 
duction; or,  C3)  that  it  allows  data  (or  analysis)  for  one  class  of 
pressure  time  histories  to  predict  the  blast-damage  results  of  other 
time  histories  (i.e.,  that  it  removes  the  problem  created  by  0 of 
equation  (6)). 

In  order  to  examine  these  questions  a number  of  different  kinds  of 
loading  histories  should  be  examined.  Physical  experiments  have  a 
limited  instrumentation  capability  and  even  greater  limitations  of 
availability.  Therefore,  it  was  necessary  to  augment  available  test 
data  with  computer  experiments  (i.e,,  structural  response  simulations). 

Six  types  of  classical  pulse  loadings  were  initially  considered  for 
this  kind  of  computer  experiment:  rectangular  (flat  top,  instant  rise 

and  fall);  quarter  cosine  (instant  rise,  cosine  fall,  no  negative  phase); 
linear  decay  (instant  rise,  linear  fall  to  zero  overpressure);  half  sine 
wave  (no  negative  phase);  ramp  (linear  rise,  instant  fall);  and,  Friedlander. 
Figure  6 shows  the  result  of  applying  equation  (2)  to  the  first  five 
types  of  loading.  This  figure  can  be  used  to  determine  the  pressure 
ratios  that  are  required  to  achieve  the  third  objective  stated  earlier 
in  this  section. 

Four  kinds  of  computer  simulations  were  used  to  obtain  data  on 
perfectly  instrumented  targets  with  which  to  study  the  proposed  model: 
single  and  five  degree-of-freedom  models;  a cylinder  model  (psuedo-Kirchoff 
theory);  and,  an  advanced  beam  vibration  model.  Figures  7 and  8 depict 
the  kinds  of  systems  simulated  by  the  first  two  models.  Figure  9 depicts 
the  time  histories  applied  to  the  one  and  five  degree-of-freedora  and  REPSIL 
models.  Figure  10  shows  the  shell  model  with  its  spacial  distribution  of 
load  having  either  rectangular  or  exponentially  decaying  time  dependent 
pressure  loads.  Figures  11  through  14  show  the  configurations,  time 
histories  and  load  parameters  of  the  advanced  beam  analysis. 


19 


20 


roposed  Model. 


RECTANGULAR  PULSE 
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Figure  9,  Pressure  Loading  Functions  Used  in  the  SWRI  and  REPSIL  Models. 


Figure  10.  Frontal  Cosine  Loading  of  Fixed-Ended  Cylinder. 


Area  Aj  _ . 9 H W;  W2 
Area  A\  ,05HWi  Wi 

Area  A3  _ .O5HW3  _ W3 
Area  Ai  .05HWi  W, 

Figure  12.  Synthesized  Beam  Cross-Section. 
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For  rh«  singit  triangU  case:  use  f2(^)  only. 

For  rhe  exponential  cose;  use  f,(r)  only 

with  n = a = 1 . 

For  the  double  triangle  case : a ^ 0 , n s i 


P,  =2Po  , = , 


Figure  13.  KADBOP  Blast  Response  Code  Forcing  Functions. 
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Figure  14. 


The  analysis  of  the  single  degree-of-freedom  system  and  the  five 
degree-of-freedom  system  were  conducted  using  a finite  difference  com- 
puter code 7 developed  by  Southwest  Research  Institute  (SWRI)  under 
contract  to  BRL.  In  that  program  the  plate  (or  beam  case  utilized) 
equations  of  motion  for  flexure  are  solved  under  assumed  moment- 
curvature  relations.  It  was  assumed  that  a bilinear  moment- curvature 
relation  with  hysteretic  recovery  could  approximate  the  response  of  a 
representative  system  made  of  aluminum  alloy.  In  Figure  8,  is  the 

elastic  rigidity,  D is  the  plastic  flexural  rigidity  and  M is  the 

bending  moment  at  which  the  transition  occurs. 

The  schematic  form  of  the  first  two  systems  investigated  are  those 
shown  in  Figure  7 where  K represents  the  spring  bending  characteristics, 

W represents  the  mass  of  each  element  and  represents  the  blast 

loading  function  to  be  applied  to  the  system.  Table  I lists  the 
parameters  used  to  analyze  those  systems. 

The  output  of  the  computer  program  was  the  maximum  bending  moment 
obtained  in  each  spring  element.  Beam  curvature  is  known  and,  assuming 
a beam  thickness,  the  resulting  maximum  strain  can  be  determined  after 
the  system  is  allowed  to  oscillate  through  enough  cycles  to  be  reasonably 
assured  that  the  maximum  strains  were  achieved.  For  the  one  degree-of- 
freedom  system  this  occurred  on  the  first  cycle  of  response.  However, 
for  the  five  degree-of-freedom  system  many  cycles  were  allowed  to  proceed 
before  stopping. 

These  two  computational  experiments  were  chosen  as  a means  of 
investigating  a wide  range  of  loading  characteristics  as  shown  in 
Figure  9. 

The  cylindrical  shell  computation  was  included  in  order  to  take 
advantage  of  the  highly  developed  REPSIL®*^  computer  code.  The  specific 
case  studied  is  a fixed  end  6061-T6  circular  cylindrical  shell  of  768 
elements  having  the  properties  listed  in  Table  I. 


7 

W.  E.  Baker y S.  SilvermcoXy  P.  A.  CoXy  Jr.  y and  D.  Young y "Structural 
Response  of  Helicopters  to  Muzzle  and  Breech  Blasty"  Vol  II y Final 
Technical  Report  SWRI  02-2029 y Contract  No.  DAAD  05-67-C-0201  (BRL) 
dated  March  1962. 

g 

N.  J.  Huffingtcriy  Jr.  y "Large  Deflection  Elastoplastic  Response 
of  Shell  Structures y " Ballistic  Research  Laboratories  Report  2575, 
November  1970 y AD  if-  717005. 

9 . 

N,  J.  Buffington y Jr. y "An  Analytical  Study  of  Explosive  Loading 
Techniques  for  Simulation  of  Impulsive  Loading  to  Structures  at 
Lethality  LevelSy  " Ballistic  Research  Laboratories  Report  No.  1621y 
November  1972y  AD  if  907436-L. 
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Table  I.  Parameters  Used  in  SWRI  and  REPSIL  Computer  Codes. 
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Length  15,24  cm  (6  in) 

Nominal  O.D.  7.62  cm  (3  in) 

Thickness  0.107  cm  (0.042  in) 


The  pressure  histories  applied  to  the  shell  act  over  the  top 
half  cylinder  having  peak  pressure  on  the  crown  line  and  a cosine 
distribution  around  to  the  half  cylinder  line  as  shown  in  Figure  10. 
The  time  histories  have  the  general  form 

P = COS0  (1  - t/t^)exp(-6t/t^)  , 


where  9 is  zero  at  the  crown  line  of  the  cylinder  and  3 = 2.0. 

The  final  system  considered  in  the  model  validation  was  a 
contractor  effort which  was  generated  under  a separate  set  of 
requirements.  It  was  included  because  it  provides  an  example  of  a 
relatively  heavy  beam  structure.  The  relatively  stiff  structure 
(clamped  at  both  ends)  filled  out  the  spectrum  of  targets  even  though 
the  loading  was  very  limited.  The  degree  of  limitation  on  the  loading 
is  made  clear  by  examination  of  Figure  14.  The  loads  are  essentially 
identical  over  60  percent  of  target  loading  regime. 

The  contractor  calculations  were  made  using  their  KADBOP^^  code. 
The  KADBOP  beam  program  employs  a finite  element  approach  wherein  the 
beam  is  represented  by  a series  of  discrete  masses  interconnected  by 
straight  weightless  bars.  Fifteen  equally  spaced  mass  points  were  used 
to  simulate  the  beam  length  and  the  beam  was  assumed  to  be  loniformly 
loaded. 


To  use  the  Kaman  data  in  the  correlated  concept  presented  here, 
it  was  necessary  to  re-compute  P”.  To  do  this  the  three  curves  of 
Figure  14  were  determined  for  the  loading  functions  used  by  the 
contractor  to  determine  non-dimensional  forms  of  the  parameters. 


cr*  peak’ 


and  F. 


VII.  CALCULATED  TEST  RESULTS 

To  test  the  results  of  the  model  we  have:  created  numerical 

simulation  of  physical  experiments;  applied  the  model  to  the  results; 
applied  the  conventional  techniques  (peak  pressure  versus  total  impulse) 
to  the  results;  and,  examined  the  relative  dispersion  of  the  two  methods. 
Visual  (i.e.,  subjective)  judgements  are  made  as  to  the  relative  merits 
of  the  two  methods.  A perfect  method  would  produce  iso-damage  curves 
that  are  independent  of  loading  pulse  shape.  Neither  system  is  perfect, 
therefore,  closeness  (minimum  distance  between)  of  curves  for  a given 
iso-damage  criterion  with  different  time  histories  is  our  standard 
of  how  good  a given  method  is.  Only  relative  measures  of  goodness 
are  possible  in  this  way. 


10 

K.  R.  WetmoTey  "An  Intewm  Teahniaal  Report  an  Phase  I of  the  BRL 
Pressure  Impulse  Blast  Program, " 30  September  1974,  Contract 
DAA0S-74-C-0742,  Kaman  Avidyne. 
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It  is  difficult  to  give  an  objective  measure  for  the  value  of  an 
engineering  practice,  but  the  authors  believe  that  the  ability  to  cause 
all  fixed  damage  iso-damage  curves  to  fall  within  experimental  scatter 
of  one  another  would  meet  the  third  objective  stated  in  Section  VI. 

In  general,  an  iso-damage  criteria  is  selected  to  satisfy  a "constant 
level -of- damage  requirement."  This  criteria  could  have  been:  maximum 

strain  anywhere  in  the  target,  maximum  deflection  of  the  target,  maximum 
deflection  at  a specified  location  on  the  target,  etc.  However,  for  the 
purposes  of  the  investigation  the  iso-damage  criteria  was  arbitrarily 
selected  to  be  the  maximum  percent  strain  on  the  surface  of  the  target  as 
calculated  at  the  discrete  grid  locations. 

Data  about  the  results  is  given  here  in  two  forms:  Tables  II  through 

V and.  Figures  15  through  24.  Tabular  data  is  included  because  the 
present  evaluation  of  the  proposed  model  is  still  subjective.  Other 
investigators  may  choose  to  use  some  tr^itional  objective  method 
of  comparing  the  data  for  both  P-I  and  P-I  iso-damage  methodology.  The 
present  authors  have  sought  the  assistance  of  experts  within  the  Ballistic 
Research  Laboratories*  to  assess  the  relative  merits  of  the  two  methods. 

No  technique  was  immediately  available;  therefore,  the  subjective 
evaluation  will  be  published  at  this  time. 

The  subjective  evaluation  used  by  the  authors  was  based  on  the 
question:  How  well  do  the  iso-damage  curves  for  one  loading  history 

(0  of  equation  (6))  predict  other  loadings?  Specifically:  How  well 

relatively,  do  the  rectangular  pulse  iso-damage  curves  (one  degree-of- 
freedom  case)  in  Figures  IS  and  16  predict  iso-damage  data  for  a linear 
decay  pulse  or  data  for  a sawtooth  pulse  or  other  type  loadings? 

Similar  questions  are  applied  to  the  five  degree-of-freedom  simulation 
data  of  Figures  17  and  18.  For  REPSIL  the  question  is  how  well  does 
rectangular  pulse  data  predict  exponential  decay  data  (Figures  19  and 
20)?  For  the  Kaman  Avidyne  beam  data  we  ask  the  reader  to  pose  his  own 
comparison  (Figures  21  through  24).  In  these  last  four  figures  the 
authors  believe  that  the  case  is  as  strongly  in  favor  of  the  P-I  method 
as  in  the  preceding  cases.  However,  it  has  been  suggested  that 
merely  displaying  a curve  fit  to  a part  of  the  data  biases  the  observers' 
opinion.  The  Kaman  Avidyne  data  was  chosen  to  be  left  without  a reference 
curve  because  it  has  the  least  variation  in  loading  function  (see 
Figures  13  and  14)  and  therefore,  subjective  judgement  of  the  results 
was  most  likely  to  be  biased  by  introduction  of  visual  reference. 

It  is  the  judgement  of  the  authors  that  visual  inspection  of 
JFigures  15  through  24  clearly  supports  the  contention  that  the  use  of 
P and  I is  superior  to  the  use  of  Ppg^j^  and  I.  The  weakest  support 

^Private  aanmuniaations : Dr.  Malcolm  Taylor 
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Table  n.  Single  Degree-of-Freedom  Iso-damage  Data 


C i-i  rsi  V 'O 


00  ^ O ltj  oi  rsj 
o ^ (N  (N 


O fsj  ^ ( o 


^ (N  (N  , 


\£>  fsj  O O 


00  CO  00  Oi  rsi  r- 


fN  fS  ^ c 


«NJ  0» 

— »«?  a.  00 
<N  ^ O K>  o 00  <J5 

fsj  im-t  o O 


tt  ^ C?i  ^ C O 
-H  -N  Oj  00 


fSJ  -N  ^ O O 


K:  hO  ^ ^ lT) 


OOCCCiO*. 


o>  o 

^ ^ r«-  C7i  Ki  CO  ^ 
TN  ••  O O Ln 


Ol  o 

— ^ h-  CT»  oe  ^ 
r^i  ^ O O 'O  LT> 


32 


Table  HI*  Five  Degree-of-Freedom  Iso-daiiage  Data 
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Table  IV.  REPSIL  Iso-damage  Data 
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Table  V 


Data  from  Kaman  Avidyne  Beam  Response  Calculations 


Strain  » 3.0%  Strain  - 15.0% 


Loading 

Average* 

Peak 

Total 

Average** 

Peak 

Total 

Function 

^essure 

Pressure 

Impulse 

Pressure 

Pressure 

Impulse 

•Type 

F(bar) 

P, 

I (bar-As) 

F (bar) 

P (bar) 

I (bar-ms) 

5.18 

0.46 

3.23 

18.61 

23.46 

11.73 

3.44 

3.99 

3.99 

11.01 

12.71 

12.71 

2.79 

3.06 

6.13 

8.61 

9.26 

18.52 

Triangular 

2.58 

2.75 

8.26 

8.09 

8.50 

25.48 

(Linear 

2.48 

2.61 

10.46 

7.87 

8.15 

32.62 

Decay) 

2.41 

2.52 

12.59 

7.74 

7.97 

39.83 

2.35 

2.42 

18.14 
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7.73 

57.98 

2.32 

2.37 

23.72 

7.49 

7.62 
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113.20 
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7.34 

367.3 

6.28 

8.60 

3.17 

23.37 

32.59 

11.99 

4.05 

5.07 

3.72 

13.29 

16.61 

12.22 

Expo- 

3.16 

3.67 

5.41 

9.86 

11.26 

16.57 

nential 

2.86 

3.21 

7.08 

8.82 

9.66 

21.33 

Decay 

2.68 

2.94 

8.67 

8.39 

9.00 

26.48 

2.59 

2.80 

10,30 

8.16 

8.62 

31.72 

2.47 

2.61 

14.39 

7.83 

8.14 

44.92 

2.40 

2.50 

18.37 

7.69 

7.92 

58.26 

2.34 

2.41 

26.57 
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6.38 

3.83 

15.34 

20.90 

12.54 

Triangular-  3.87 

4.63 

5.56 

12.40 

14.74 

17.69 

Triangular 

3.39 

3.95 

7.10 

10.71 

12.36 

22.24 

(Bilinear 

3.10 

3.52 

8.44 

9.61 

10.73 

25.75 

Decay) 

2.92 

3.23 

9.75 

8.96 

9.79 

29.37 
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2.88 
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16,19 
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* Pj.j.  = 2.21  bar 

**  Per  =»  7.24  bar 
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Figure  16.  Plot 
One 


•3.83 


2 4 6 8 10  12  14  16  18 

Tofal  Impulx  (Bar -ms) 


Figure  17.  Plot  of  Vs.  Total  Impulse  for 

Five  Degree-of-Freedom  Data. 


38 


VI  I I I 1 1 1 1 1 I 

0 2 4 6 6 10  12  14  16  18 

Total  Impulsa  (8ar-ms) 


Figure  18.  Plot  of  P Vs.  Total  Impulse  for 
Five  Degree-of-Freedom  Data. 
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Figure  19.  Plot  of 
REPSIL  1 


I 


Vs.  Total  Impulse  for 
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Figure  22,  Plot  of  P Vs.  Total  Impulse  For  Kaman  3 Percent  Strain  Data 


o 
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Figure  23.  Plot  of  Ppg^l^  Vs.  Total  Impulse  for  Karaan  15  Percent  Strain  Data. 
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Figure  24,  Plot  of  P Vs.  Total  Impulse  for  Kaman  15  Percent  Strain  Data 


for  this  position  comes  from  the  contractor  data.  This  is  a consequence 
of  the  strong  similarity  of  the  pulse  shapes  used  (recall  that  all  three 
are  approximations  to  the  same  time  history).  No  case  shows  the 
traditional  method  to  be  superior  to  the  model  proposed  here. 


VIII.  SHOCK  PULSE  EFFECTIVENESS 

In  making  comparisons  of  the  pulses  parameterized  in  Figures  6 and 
14  one  may  ask:  How  are  the  traditional  model 

^*^peak  “ *^cr^^^total  ~ ^cr^  “ ^t 

and  the  proposed  model 

(P  - P ^){I  ^ - I ) = K 

cr-^  ^ total  cr-^ 

related?  Looking  at  Figure  6 one  sees  that  ^Ppg^  is  almost  a linear 
function  of  Pgj./Ppgak*  Making  use  of  this  characteristic  to 
approximate  "P/P  ^ by 

Uw  3JC 


peak 

we  find  that 


= + (1  - K^)  p 


cr 


peak 


P = K.P  , + (1  - KJP 

1 peak  F cr 


which  may  be  substituted  into  the  new  iso-damage  model  with  the  result 
* ■’cr  - ■=l''cr  ‘ ’’cr’ ("total  ' "cr"  = 

This  simplifies  to  give 


(P. 


- P„J(I. 


- I J = 


peak  ‘ cr"" '■"total  ‘cr''  K 


Thus  in  the  linear  approximation  to  P/P  , we  find  that  the  P /P 

peak  cr'  peak 

intercept,  K^,  is  the  transformation  factor  between  the  traditional 
and  new  models,  i.e.,  = K/Kj^) . In  other  words 
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•=!  * ['Vak  - ('total  - 'cT»  = C ' ’’or’ "total  ‘ 'or’  ' >= 

Thus,  we  see  that  the  proposed  model  is,  in  effect,  a means  of  modifying 
the  material  contained  in  the  traditional  P-I  maps.  In  general,  it  can 
be  said  that  the  parameter,  is  a quantitative  measure  of  the 

effectiveness  of  a pulse  as  compared  to  a rectangular  pulse  having  the 
same  total  impulse  in  causing  damage  irrespective  of  the  target 
characteristics. 


IX.  USE  OF  TOE  PROPOSED  MODEL  IN  DATA  REDUCTION 


So  far,  the  discussion  has  been  limited  to  the  theoretical  advan- 
tages of  the  present  model.  Another  aspect  which  is  of  equal  if  not 
greater  importance  is  that  the  proposed  model  is  a practical  tool  for 
the  reduction  of  test  data.  In  short,  the  discussion  so  far  has 
neglected  objectives  one  and  two  of  Section  VI. 


Figure  25  shows  what  we  can  expect  to  gather  in  the  way  of  data 
from  blast  test  instrumentation.  As  one  can  see,  the  loading  is  not 
classical.  In  applying  methods  that  depend  on  peak  pressures,  total 
impulse  and  similarity  to  classical  loading,  the  data  reduction  relies 
heavily  on  subjective  judgement.  With  the  proposed  model,  the  results 
are  integrals  of  time  histories  and  as  such  are  not  dependent  on  sub- 
jective evaluation  of  empirical  records.  It  will  be  possible  to  use 
the  ^“I^otal  pair  or  its  derivative  characterizations  of  the 

target  response  to  a given  attack  in  either  theoretical  or  experimental 
studies. 


Use  of  this  model  for  data  analysis  should  incur  no  additional  cost 
because  most  agencies  recording  tests  data  now  have  the  pressure-time 
histories  digitized  for  computer  plotting  and  impulse  determination. 
Figure  25  provides  a case  in  point.  Data  for  two  of  five  side-on  pres- 
sure measurements  taken  on  one  Fuel  Air  Explosive  (FAE)  test  are  shown. 
The  FAE  device  was  initiated  near  the  ground  and  identified  as  having  a 
nominal  charge  weight  of  33,6  kg. 

In  Figure  26  we  see  the  damage  number  CDN)  versus  radius  for 
this  event.  Also  included  on  the  figure  are  the  damage  numbers  versus 
radius  for  bare  spherical  pentolite  charges  detonated  in  close  proximity 
to  the  ground  and  assumed  to  have  an  effective  charge  weight  of  1.8 
times  actual.  The  conclusion  that  we  draw  from  this  comparison  is 
that  for  side-on  pressure  loading  cases  outside  the  cloud  and  near 
the  ground,  this  FAE  device  is  approximately  2.1  times  as  effective  as 
Pentolite.  This  indicates  that  relatively  soft  targets  can  have  damage 
numbers  up  to  a limiting  value  as  indicated  by  the  knee  in  the  curve. 
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OVERPRESSURE  (bar)  OVERPRESSURE  (bar) 


33.6Kq  ETHYLENE  OXIDE 


PRESSURE  AT  6.1  METRES 


33.6 Kq  ETHYLENE  OXIDE  PRESSURE  AT  10.3  METRES 


Figure  25.  Typical  Fuel-Air  Explosive  Pressure  Time  Histories. 
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DN  (bar^  sec) 


RADIUS  (METRES) 


Figure  26.  Damage  Number  Vs.  Radius  for  Near  Surface 
Bursts  of  Pentolite  and  Fuel-Air  Explosive. 
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Another  example  of  the  use  of  this  model  in  interpreting  test 
data  is  shown  in  Figure  27.  This  figure  shows  localized  loading  time 
histories  for  two  helicopters  exposed  to  a large-scale  high-explosive 
detonation. The  difference  between  the  two  curves  results  from  one 
helicopter  being  revetted  (lower  curve)  and  the  other  one  being  un- 
revetted. One  can  see  that  the  determination  of  the  damaging  potential 
of  the  loading  functions  (which,  incidentially , for  this  case  is  a 
measure  of  the  ability  of  the  revetment  to  reduce  damage)  is  enhanced 
as  judged  by  use  of  the  damage  number  criterion.  Damage  numbers  for 
these  two  loading  functions  are  included  on  the  plots.  It  can  be  seen 
that  the  damage  momber  for  the  non-revetted  helicopter  is  slightly 
larger  than  that  for  the  revetted  helicopter.  The  ratios  of  the 
damage  numbers  (HELO  2/HELO  3)  was  found  to  be  1,14. 

If  one  were  to  use  the  previous  concepts  for  evaluating  these  same 
data,  one  would  get  a peak  pressure  ratio  of  1.19,  a total  impulse 
ratio  of  1.05  and  a peak  pressure-impulse  product  ratio  of  1,26.  These 
numbers  indicate  that  the  results  would  depend  upon  the  ratio  or  ratios 
considered  important.  It  is  interesting  to  note  that  the  average  of 
these  three  ratios  is  1.16,  which  is  close  to  the  ratio  obtained  using 
our  model.  It  thus  appears  that  our  model  can  be  of  significant  help 
in  evaluating  test  results  as  well  as  in  quantifying  test  data. 


X,  CONCLUSIONS 

A conclusion  which  has  been  a defacto  part  of  the  blast  community 
is  that  P-I  methods  are  not  applicable  in  the  pure  impulse  and  static 
load  regions  of  dynamic  structural  response  analysis.  This  conclusion 
has  been  brought  into  sharp  focus  by  comparison  of  the  presently  pro- 
posed P-I  method  with  its  traditional  and  less  versatile  predecessor. 

The  proposed  method  provides  an  explanation  for  interesting  previous 
methods  having  differing  degrees-of-sophistication.  These  logical 
connections  relate  the  one  parameter  method  of  0.  T.  Johnson  to  the 
two  parameter  P-I  method.  In  turn,  the  P-I  method  is  related  to  experi- 
ment and  theory  in  structural  blast  damage. 

The  currently  proposed  method  provides  a convenient  standard  method 
of  blast  data  reduction  which  may  be  applied  to  future  experimental 
work  and  to  correlate  future  data  to  that  existing  today. 

The  proposed  method  provides  a means  of  relating  damage  from  explo- 
sives having  known  characteristics  to  that  caused  by  explosives  having 
significantly  different  blast  characteristics. 

Finally,  the  present  research  has  lead  to  an  ongoing  effort  to 
apply  statistical  methods  to  validation  of  P-I  iso-damage  methodology. 
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OVERPRESSURE  (bof)  OVERPRESSURE  (bor) 


MEASUREMENT  NO.  P2  AT  426  METRES 
CHARGE  WT*0,3KT  TNT 


HELICOPTER  NO.  2 


Figure  27.  Pressure  Loading  Functions  on  Two  Helicopters 
During  a Large  Scale  Blast  Field  Trial. 
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1 AVCO  Corporation 

Structures  and  Mechanics  Dept. 
ATTN:  Dr.  William  Eroding 

Mr.  J.  Gilmore 
Wilmington,  MA  01887 

1 The  BDM  Corporation 

ATTN:  Technical  Library 

1920  Aline  Avenue 
Vienna,  VA  22180 


No.  of 

Copies  Organization 

1 Bell  Telephone  Labs,  Inc, 

ATTN:  Tech  Rpt  Ctr 

Mountain  Avenue 
Murray  Hill,  NJ  07971 

1 The  Boeing  Company 

ATTN : Aerospace  Library 

P.O.  Box  3707 
Seattle,  WA  98124 

2 California  Research  and 
Technology,  Inc. 

ATTN:  Ken  Kreyenhagen 

Technical  Library 
6269  Variel  Avenue 
Woodland  Hills,  CA  91364 

1 Calspan  Corporation 

ATTN:  Technical  Library 

P.O.  Box  235 
Buffalo,  NY  14221 

1 Civil/Nuclear  Systems 
Corporation 

ATTN : Robert  Crawford 

1200  University  N.E. 
Albuquerque,  NM  87102 

1 EG§G,  Incorporated 
Albuquerque  Division 
AITN:  Technical  Library 

P.O.  Box  10218 
Albuquerque,  NM  87114 

1 General  American  Trans 
Corporation 

General  American  Research 
Division 

ATTN:  G.  L.  Neidhardt 

7449  N.  Matchez  Avenue 
Niles,  IL  60648 

1 General  Electric  Company-TEMPO 
ATTN:  DAS I AC 

P.O.  Drawer  QQ 
Santa  Barbara,  CA  93102 
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I President 

General  Research  Corporation 
ATTN : Library 

McLean,  VA  22101 

1 J.  G.  Engineering  Research 
Associates 
3831  Menlo  Drive 
Baltimore,  MD  21215 

1 J.  L.  Merritt 
Consulting  § Special 

Engineering  Services,  Inc. 
ATTN:  Technical  Library 

P.O.  Box  1206 
Redlands,  CA  92373 

2 Kaman  AviDyne 

ATTN:  Dr.  N.  P.  Hobbs 

Mr.  S.  Criscione 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01803 

2 Kaman  Sciences  Corporation 
ATTN : Library 

Dr.  Don  Sachs 

1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80907 

2 Lockheed  Missiles  5 Space  Co. 
ATTN:  T.  Gerrs  D/52-33, 

Bldg.  205 

Tech  Info  Ctr, Doc/ Coll . 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

1 McDonnell  Douglas  Corporation 
ATTN:  Robert  W.  Halprin 

5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

1 Honeywell,  Inc. 

Govt  5 Aero  Products  Div 
600  Second  Street  N.E. 

Hopkins,  MN  55343 


No . of 

Copies  Organization 

1 The  Mitre  Corporation 

ATTN : Library 

RT  62  and  Middlesex  Turnpike 
P.O.  Box  208 
Bedford,  MA  01730 

2 Physics  International  Corp. 

ATTN:  E.  T.  Moore 

Dennis  Orphal 
2700  Merced  Street 
San  Leandro,  CA  94577 

4 Physics  International  Corp. 

ATTN:  Robert  Swift 

Charles  Godfrey 
Larry  A.  Behrmann 
Technical  Library 
2700  Merced  Street 
San  Leandro,  CA  94577 

5 RSD  Associates 

ATTN:  Dr.  H.  L.  Brode 

Dr.  Albert  L.  Latter 
Bruce  Hartenbaum 
William  B.  Wright 
Henry  Cooper 
P.O.  Box  9695 
Marina  del  Rey,  CA  90291 

4 R^D  Associates 

ATTN:  Jerry  Carpenter 

Sheldon  Schuster 
J.  G.  Lewis 
Technical  Library 
P.O.  Box  9695 
Marina  del  Rey,  CA  90291 

1 Sandia  Laboratories 

ATTN:  Doc  Control  for  3141 

Sandia  Report  Coll. 
P.O.  Box  5800 
Albuquerque,  NM  87115 

1 John  Hindes,  Consultant 
P.O.  Box  6434 
Santa  Barbara,  CA  93111 
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Organization 


2 Sandia  Laboratories  2 Terra  Tek,  Inc. 

Livermore  Laboratory  ATTN:  Sidney  Green 

ATTN:  Doc  Control  for  Tech  Lib  Technical  Library 


Doc  Control  for  L.  Hill 
P.O.  Box  969 
Livermore,  CA  94550 

1 Electromechanical  Systems  of 

New  Mexico,  inc. 

ATTN:  R.A.Shunk 
PO  Box  11750 
Albuquerque,  NM  87112 

3 Science  Applications,  Inc. 
ATTN:  R.  Seebaugh 

William  M.  Layson 
John  Mansfield 
1651  Old  Meadow  Road 
McLean,  VA  22101 

2 Science  Applications,  Inc. 

ATTN:  David  Bernstein 

D.  E.  Maxwell 
7850  Edgewater  Drive 
Oakland,  CA  94621 

2 Science  Applications,  Inc. 
ATTN:  Technical  Library 

Michael  McKay 
P.O.  Box  2351 
LaJolla,  CA  92038 

1 Shock  Hydrodynamics,  Inc. 

A Division  of  Whittaker  Corp. 
ATTN : L . Zernow 

4710-16  Vineland  Avenue 
North  Hollywood,  CA  91602 

5 Systems,  Science  § Software 
ATTN;  Robert  T.  Allen 
Donald  R.  Grine 
Ted  Cherry 
Thomas  D.  Riney 
Technical  Library 
P.O.  Box  1620 
La  Jolla,  CA  92037 


420  Wakara  Way 

Salt  Lake  City,  UT  84108 

2 Tetra  Tech,  Inc. 

ATTN:  Li -San  Hwang 

Technical  Library 
630  North  Rosemead  Blvd. 
Pasadena,  CA  91107 

7 TRW  Systems  Group 

ATTN:  Paul  Lieberman 

Benjamin  Sussholtz 
Norm  Lipner 
William  Rowan 
Jack  Farrell 
Pravin  Bhutta 
Tech  Info  Ctr/S-1930 
One  Space  Park 
Redondo  Beach,  CA  90278 

2 TRW  Systems  Group 

ATTN:  Mr.  F.  A.  Pieper 

Greg  Hulcher 

San  Bernardino  Operations 

P.O.  Box  1310 

San  Bernardino,  CA  92402 

2 Union  Carbide  Corporation 

Holifield  National  Laboratory 
ATTN : Doc  Control  for 

Tech  Lib 

Divil  Defense  Research 
Proj  . 

P.O.  Box  X 

Oak  Ridge,  TN  37830 

1 Universal  Analytics,  Inc. 
ATTN:  E.  I.  Field 

7740  W.  Manchester  Blvd. 

Playa  Del  Rey,  CA  90291 
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1 Massachusetts  Institute  of 

Technology 

Aeroelastic  and  Structures 
Research  Laboratory 
ATTN:  Dr.  E.  A.  Witmer 

Cambridge,  MA  02139 

2 Southwest  Research  Institute 

ATTN : Dr . W . E . Baker 

A.  B.  Wenzel 
8500  Culebra  Road 
San  Antonio,  TX  78206 

4 Stanford  Research  Institute 
ATTN:  Dr.  G.  R.  Abrahamson 

Carl  Peterson 
Burt  R.  Gasten 
SRI  Library,  Rin  G02I 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

I University  of  California 
Berkeley  Campus,  Rm  8 
ATTN : G . Sackman 

2543  Channing  Way 
Berkeley,  CA  94720 

1 University  of  Dayton 
ATTN:  Hallock  F.  Swift 

300  College  Park  Avenue 
Dayton,  OH  45409 

1 University  of  Illinois 

Consulting  Engineering  Services 
ATTN:  Nathan  M.  Newmark 

1114  Civil  Engineering  Building 
Urgana,  IL  61801 

2 The  University  of  New  Mexico 


No . of 

Copies  Organization 

1 URS  Research  Company 

ATTN:  Technical  Library 

155  Bovet  Road 

San  Mateo,  CA  94402 

1 Weidlinger  Assoc.  Consulting 
Engineers 

ATTN:  J.  Isenbert 

2710  San  Hill  Road 
Suite  230 

Menlo  Park,  CA  99025 

1 Westinghouse  Electric  Company 
Marine  Division 
ATTN:  W.  A.  Volz 

Hendy  Avenue 
Sunnyvale,  CA  94008 

1 Battelle  Memorial  Institute 

ATTO : Technical  Library 

505  King  Avenue 
Columbus,  OH  43201 

4 Denver  Research  Institute 
University  of  Denver 
ATTN:  Mr.  J.  Wisotski 

Fred  P.  Vanditti 
Ron  W.  Buchanon 
Technical  Library 
P.O.  Box  10127 
Denver,  CO  80210 

3  IIT  Research  Institute 
ATTN:  Milton  R.  Johnson 

R.  E.  Welch 
Technical  Library 
10  West  35th  Street 
Chicago,  IL  60616 

2 Lovelace  Foundation  for 

Medical  Education 
ATTN:  Asst  Dir  of  Research/ 

Robert  K.  Jones 
Technical  Library 
5200  Gibson  Blvd,  SE 
Albuquerque,  NM  87108 


The  Eric  H.  Wang  Civil 

Engineering  Research  Facility 
ATTN:  Larry  Bickle 

Neal  Baum 
University  Station 
Box  188 

Albuquerque,  NM  87131 
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2 Washington  State  University 
Administration  Office 
ATTN:  Arthur  Miles  Hohorf 

George  Duval 
Pullman,  WA  99163 


Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Dir,  USAMSAA 

ATTN:  Dr.  J.  Sperrazza 

Mr.  R.  Norman,  GWD 
Dr.  Rivello 
R.  Bailey 
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